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I show that entanglement between two qubits can be generated if the two qubits interact with a 
common heat bath in thermal equilibrium, but do not interact directly with each other. In most 
situations the entanglement is created for a very short time after the interaction with the heat bath 
is switched on, but depending on system, coupling, and heat bath, the entanglement may persist 
for arbitrarily long times. This mechanism sheds new light on the creation of entanglement. A 
particular example of two quantum dots in a closed cavity is discussed, where the heat bath is given 
by the blackbody radiation. 



Since the discovery of quantum mechanics, "entangle- 
ment" has been considered a hallmark of quantum behav- 
ior [1] . Two quantum systems A and B in a pure state are 
called entangled, if their quantum mechanical state vec- 
tor can not be written as product of two states \4>a) 
and \4>b) in the Hilbert spaces of A and B, respectively. 
The last few years of research on quantum information 
processing have lead to the picture of entanglement as a 
precious resource. Entanglement plays an important role 
in super-dense coding [2] and quantum teleportation [3] , 
and is necessary for the exponential speed-up of quantum 
algorithms compared to classical algorithms [4]. 

Recently investigated examples of the controlled cre- 
ation of entanglement include trapped ions that interact 
electrostatically (or more precisely exchange phonons in 
a chain of ions [5]), and the entanglement of atoms in a 
cavity by the interaction with a specific electromagnetic 
mode of the cavity [6,7]. In the latter example entan- 
glement can be created even in the case where the cavity 
mode is itself coupled to many more degrees of freedom of 
the electromagnetic environment, i.e. if the cavity is more 
or less leaky. Nevertheless, in all these examples a third 
system with one or few degrees of freedom is clearly sin- 
gled out by mediating the interaction between the atoms 
or ions. This is true even for strongly leaking cavities 
capable of supporting super-radiance [8] , which may still 
entangle atoms [9]. 

In the following I show that entanglement can be cre- 
ated if the two systems interact neither directly, nor with 
a third system with only one or a few singled out de- 
grees of freedom, but interact with the (possibly infinitely 
many) degrees of freedom of a heat bath in thermal equi- 
librium. This is a priori not obvious since interactions 
with a heat bath lead typically to very rapid decoher- 
ence [10], thus to classical states and the destruction of 
quantum entanglement. Indeed, we will see that the en- 
tanglement created may die again on a decoherence time 
scale of the system. However, notable exceptions exist: 
i) if the two systems are coupled in a symmetric way 
to the environment the entanglement will be protected 
[11] — much in the spirit of what is known from coher- 



ent rotational tunneling [12], long living Schrodinger cat 
states [13] and decoherence free subspaces [7,14,15]. ii) 
Many environments will lead, for systems with degener- 
ate energy levels and in sufficiently high dimension, to in- 
complete decoherence, a surprising effect to be discussed 
below. 

When dealing with a "heat bath" , i.e. another system 
with very many degrees of freedom over which we do not 
have microscopic control, the definition of entanglement 
has to be generalized to mixed states. A state of a bi- 
partite system is said to be "separable", iff the density 
matrix of the state can be written as 
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where the pi are probabilities (0 < pi < 1), pf and pf are 
density matrices for the two subsystems A and B, and N 
is an arbitrary integer. A state that is not separable is 
said to be entangled [16]. A simple criterion for entangle- 
ment of bipartite systems of dimensions 2 x 2 or 2 x 3 was 
proven by the Horodecki family [17]: A state p of a 2 x 2 
or 2 x 3 bipartite system is separable, iff p has a non- 
negative partial transpose p TA . The partial transpose 
p TA is obtained by transposing in a matrix representa- 
tion of p only the indices corresponding to subsystem A, 

Le - Pik,jl = Pjk,il with p 3 k,u = (j| <S> (Mp\i) <S> 10- 

Suppose Alice and Bob both own a qubit with ba- 
sis states |0) and |1) over which they have local con- 
trol. The qubits do not interact with each other. Thus, 
the Hamiltonian representing the two qubits is simply 
Hab = H\ + Hb , where Ha acts only on Alice' Hilbert 
space, and Hb only on Bob's. Suppose further that the 
qubit states |0) and |1) are energy eigenstates with de- 
generate energies, for both qubits. In this case, Hab = 
up to an irrelevant constant. For situations where the 
degeneracy is not exact, let us assume at least that the 
inverse level spacing is much larger than any time scale 
that we are interested in. The dynamics induced by a 
finite Ha or Hb can then be neglected and we can drop 
the "system Hamiltonian" Hab [18]. The heat bath will 
be described as a collection of N harmonic oscillators, 
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For the interaction with the heat bath we assume a cou- 
pling Hamiltonian 



Hint = (S A + S B )B , B^J29kQk- 



(3) 



where S A and S B are "coupling agents" acting on the 
Hilbert spaces of Alice and Bob, respectively, and the gu 
are coupling constants to the fcth oscillator. An example 
that is described by (3) will be analyzed in detail below. 

Let us further assume that the qubit basis states |0) 
and |1) are eigenstates of S A and S B with eigenvalues 
ao, Oi, and bo, &i, respectively. The combined compu- 
tational basis states |00), |01), |10), and |11) are then 
eigenvectors of S A + S B with corresponding eigenvalues 
o-o + bo, «o + bi, a\ + bo, and a\ + &i, respectively. 

Protecting their qubits momentarily from the environ- 
ment, Alice and Bob prepare pure initial states \<p A ) and 
\p B ) of their respective qubits. The heat bath is assumed 
to be initially in thermal equilibrium at temperature T, 
and so the total initial state is the density matrix 



W(0) = \<p A )(tp A \®\tp B )(tp B \®L 



-ffbath/fcBT 
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where Z = trb a the~ ffbath ^ fcBT and fee denotes Boltz- 
mann's constant. The time evolution of Alice' and Bob's 
qubits alone is described by the reduced density matrix 
p[t) = tib&thW{t). That time dependence was calculated 
for an arbitrary system with negligible system Hamilto- 
nian and coupled as in (3) to a heat bath of harmonic 
oscillators in [18]. The result can be phrased in terms of 
two functions f(t) and <p(t), 
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where denotes the thermal occupation of the fcth mode 
and C(t) = (B(t)B(0)) represents the thermal bath cor- 
relation function. In the basis of eigenstates of S A + S B 
(the "pointer basis" [10]), the time evolution of pij.ki(t) 
{i,j,k,l = 0,1) is given by 

Pij,kl(t) = Pij,kl(°) ex P (-(«» + b i - Ofc - b i) 2 f(t) (7) 
+i((ai + bjf - (a fc + bif)^{t) 



In general this time evolution leads to a rapid decay of 
the off-diagonal matrix elements — unless S A + S B has 
degenerate eigenvalues. 

Suppose Alice and Bob prepare the initial states 
W A ) - (|0) - |1))/V2 and \cf> B ) = (|0) + \l))/y/2, i.e. 



p(0) = 
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and assume for the moment the symmetric coupling sit- 
uation ao = bo = 0, and a\ = b\ = 1, absorbing eventual 
prefactors into the coupling constants gk- It is straight 
forward to compute numerically the eigenvalues of the 
partially transposed density matrix p TA for given f(t) 
and (p(t). Note that these functions vanish at t = 
and are strictly positive for times t > 0; for small t 
(u>kt <C 1) always both f(t) and g(t) become finite, with 
f(t) cx ip(t) 2 / 3 . By parameterizing the eigenvalues di- 
rectly by f(t) and <p(t) one can examine all possible (har- 
monic) heat baths at the same time. A given heat bath 
leads to a certain path in the /, ip plane. Fig.l shows the 
smallest eigenvalue A of p TA as function of f(t),ip(t). 
The eigenvalue is zero for t = 0, where both / and tp 
vanish: since the two qubits were prepared in a product 
state, the partial transpose is the same as the original 
matrix, and the Schmidt decomposition gives one eigen- 
value unity and three equal zero. As soon as f(t) and 
ip(t) aquire a finite value, Ao becomes negative, however, 
meaning that the two qubits get entangled. For larger 
values of / and tp, the absolute value of Ao decays again, 
and asymptotically, for f(t) — ► oo, the state 



(9) 



is reached, independent of the behavior of the imaginary 
part. Alice's partial transpose p^£ of this matrix has 
eigenvalues 1/2, 1/4 (doubly degenerate) and zero, so 
that for f(t) — > oo no entanglement is left. Using sec- 
ond order perturbation theory in the deviation of p{t) TA 
from p^£ one easily shows, however, that for all arbitrar- 
ily large but finite f(t) the two qubits stay entangled. 

Finite entanglement is created at short times also in 
the case that ao + b\ and a\ + bo arc not degenerate. 
Numerical investigation shows that for a given /(<) and 
tp{t) the positivity of p TA may be even more strongly 
violated for non-perfect degeneracy. Non-perfect degen- 
eracy changes things drastically, however, for large f(t) 
when for all finite deviations from degeneracy the non- 
entangled state poo = jdiag(l, 1, 1, 1) is reached. There- 
fore there will be a finite time after which the initial state 
becomes separable — if f(t) reaches large values. 
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It is easy to show numerically that the choice of the 
initial state is not crucial. As long as both states con- 
tain components of both |0) and |1), the heat bath 
creates entanglement between the two qubits. And 
I have also checked that the interaction with a com- 
mon heat bath can create entanglement between a 
qubit and a qutrit (i.e. a 2 x 3 bipartite system). 



imation. This is suitable for temperatures where only 
modes with wave lengths much larger than d are popu- 
lated, i.e. ksT <C 2nhco/d (cq is the speed of light). The 
interaction in dipole approximation reads 



jjA.B _ 
H int ~ 



V P?J B \i)(j\^e k (x AtB )q k , (10) 



D A,B 
i,j=0 




FIG. 1. Smallest eigenvalue of p TA as a function of / and <p 

Let me finally propose a concrete system where the 
effect might in principle be observable. Consider two 
double-well quantum dots enclosed in an ideally conduct- 
ing, box-shaped cavity, with edge dimensions a, b, and c 
in x, y, and z directions, respectively. The two quantum 
dots are assumed identical, with the two-dimensional 
electron gas in the y — b/2 plane, and with two iden- 
tical wells to the right and left (in z direction) each, 
which might be electrostatically defined by suitable gate- 
electrodes. The symmetry centers of the double-well 
quantum dots are located at xa = (a/4, 6/2, c/2) and 
xb = (3o/4, b/2, c/2) for dot A and B, respectively. The 
centers of the wells are separated by a distance d and 
we will assume that there exist in both dots two states 
|0) and |1), localized in the right and left well, such that 
they are eigenstates of the dipole operators, i.e. er^, 
(0|r A |0) = -d/2(0,0,l) = -(l|r A |l), (0|r A |l) = for 
dot A, where e is the electron charge and ya = x — x^ 
is the position of an electron with respect to the center 
of the dot. This can be achieved to very good approx- 
imation by a very high barrier between the two wells, 
which leads to exponentially small overlap of the two 
wave functions, and negligible tunneling splitting. For 
dot B the two states are chosen in the opposite way, 
(0|r B |0) = d/2(0,0,l) = -(l|r B |l). The cavity sup- 
ports TE and TM modes ("transversal" relative to the 
arbitrarily chosen z-direction as propagation direction). 
For the above geometry the two dots interact only with 
the TM modes, if we describe the interaction between 
the dots and the electromagnetic field in dipole approx- 



where Pfj B = a,b\]) are the dipole matrix elements 
defined above, £k = ^J^j^i'fix, <Pyi fz) (with eq and /ito 
electric permeability and magnetic susceptibility of vac- 
uum, in SI units), and is the electric field amplitude 
of mode k with the dimension of a length, chosen as co- 
ordinate of the harmonic oscillator in the quantization of 
the field [19]. The mass m introduced formally for this 
purpose will cancel out again in the final expressions for 
/(f) and tp(t). The functions <p x , ip y , and ip z define the 
spatial structure of the modes. Here, only tp z is needed, 
which for perfectly conducting walls is given by [20] 
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sin(fc x x) sin(fcj,y) cos(fc z z) ; (11) 
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V = abc denotes the total volume, k± = ^Jk 2 

are the transverse wave numbers, and the wave vector 
is given by {k Xl k y ,k z ) = ir(m/ 'a,n/b,p/ 'c), m,n,p — 
0,1,2,.... Thus, 



rjA,B _,_ A.B A,B - 

H int =T<V 2^ 9k 
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/'k = : ' l \l-j^^shi{k x ^)sm(k y ^)cos(k z ^) (13) 



B , I m kx . ,, 3a, . ,, b. . c, . 
ffk = ed \ — F7 — sm(k x — )sm(k y -)cos(k z -). (14) 

The upper sign refers to dot A, the lower to B, and 
the operators a x ,B are defined for both systems as 
&x = |0)(0| — |1)(1|- One easily sees that the modes 
with even m couple to a x + a x , those with odd m cou- 
ple to a x — erf . However, the odd modes can be sup- 
pressed by a very thin, uncharged, and perfectly conduct- 
ing wire in the z direction along the x — a/2,y = b/2 
axis of the cavity, since they have non-vanishing tan- 
gential electrical field at the position of the wire. We 
then obtain the coupling Hamiltonian (3) with S A = <r x , 



S B = at 



,9k 



— g£, and quantized wave vectors k 



7r((4n x + 2)/a, (2n y + l)/b, 2n z /c), n x ,n y ,n z = 0, 1,2, .. .. 
The resulting expressions for /(f) and ip(t) are divergent, 
and the sum over k needs a cut-off. For a cavity made out 
of a real metal a natural cut-off frequency is the plasma 
frequency u> p , since the metal looses its reflectivity for 
lj > Lu p [21]. Converting the sums over k into integrals, 
defining ((d,0) = e 2 d 2 ^ / '(tt 2 c h 3 (3 2 ) , and r = ph/2, 
one finds /(f) = ((d,/3)f(t) and ip(t) = ((d,l3)<p(t) with 
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where c{x) is a cut-off function (to be specific, say c(x) — 
exp(— x/xmax)) For an aluminum cavity, htu p = 15.3eV 
[21], and we have, at T = lOOmK, r ~ 3.810" n s and 
i mEU[ = ujpT ~ 8.8 10 5 . The main contribution to the 
integrals therefore stems from x ^> 1, where we can ap- 
proximate the coth-function by 1. For the exponential 
cut-off we then have 

cos (2 arctan(te max /r)) ' 
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a function that saturates for fx max ^ r at the value x^^, 
after reaching a maximum at tjr ~ a:"^ . For quantum 
dots with d = lOnm, T = lOOmK, ((d,f3) ~ 1.8 10" 15 , 
and f(t) reaches a maximum after a time of the or- 
der 10~ 17 s before saturating at /(f) ~ 0.0014. This 
means that decoherence remains incomplete even for 
non-symmetric couplings, and the entanglement created 
by the interaction with the heat bath is preserved, till 
other decoherence mechanisms neglected in the above 
analysis kick in. Note that such incomplete decoherence 
is a rather general result for systems with degenerate en- 
ergy levels. In fact, by integrating the time dependent 
part in eq.(15) from zero to t one obtains for t — > oo a 
Dirac delta function at x = 0, and the remaining integral 
over x will give a finite constant. Thus, the time depen- 
dent part of f(t) has to decay for t — ► oo faster than 1/t, 
leaving the time independent part J dx x coth(x)c(x) / A. 
Note that the factor x, the spectral weight of the heat 
bath at zero frequency, is essential in this reasoning. De- 
coherence will always remain incomplete (in the sense 
of finite /(f) for f — > oo, depending on the circumstances 
even /(f) <C 1) between degenerate energy levels for spec- 
tral weights that vanish at zero frequency faster than the 
first power of the frequency. 

The presented scheme has an advantage over conven- 
tional creation of entanglement if Alice and Bob are so 
far apart that a direct interaction is difficult to achieve. 
Since /(f) and tp(t) do not depend on the volume of the 
cavity, very large cavities should be possible with corre- 
sponding large distances between Alice and Bob; using 
the non-symmetric coupling scheme one might even en- 
visage to dispose of the cavity altogether and rely on 
the long wave-length continuum of the cosmic electro- 
magnetic background radiation to create entanglement 
between very remote quantum dots. While only a small 
amount of mixed state entanglement will be created, it 
is well known that all entanglement of a 2 x 2 bipartite 
system can be distilled into a pure entangled state [22], 
given sufficiently many realizations of the input states 
and local coherent control. 

As a conclusion I have shown that entanglement can 
be created between two qubits that interact solely with 



a common heat bath with very many degrees of freedom. 
The explicit example of two double-well quantum dots in 
a cavity was calculated, and the phenomena of "incom- 
plete decoherence" was revealed, which may, as much 
as symmetric couplings to the environment, preserve the 
entanglement created by the heat bath. 
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